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Abstract. While the fluorescence and the ground counter techniques for the detection of ultra-high energy cosmic rays 
(UHECR) were being developed for decades, the interest in the radio detection diminished after the initial experiments in 
the 1960s. As a result, the fluorescence and the surface array techniques are more mature today, providing more reliable 
measurements of the primary cosmic particle energy, chemical composition and the inelastic cross-section. The advantages of 
the radio technique are 100% duty cycle and lower deployment and operational costs. Thus, the radio technique can greatly 
complement the fluorescence and the ground array detection and can also work independently. With the ANITA balloon 
detector observing UHECRs and the success of LOPES, CODALEMA and other surface radio detectors, the radio technique 
received a significant boost in recent years. Reliable Monte Carlo (MC) simulations are needed in order to obtain the energy 
and other parameters of the primary cosmic ray particle from the radio observations. Several MC techniques, like ZHairesS 
and the Endpoint Formalism, were proposed in recent years. While they seem to reproduce some of the observed data quite 
well, there is a divergence between the different approaches under certain conditions. In this work we derive these approaches 
from Maxwell's equations and prove their identity under certain conditions as well as discuss their applicability to the UHECR 
air showers and to a proposed experiment at SLAC. 
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I. INTRODUCTION 

Ultra high energy cosmic rays (UHECRs) are charged 
particles with the energy above 10 18 that are constantly 
bombarding the earth. Even after decades of intensive 
study, the origin, acceleration mechanism and chemical 
composition of the UHECRs are not yet known. While 
a correlation between the UHECR arrival directions and 
the nearby active galactic nuclei has been reported [1], 
the correlation becomes weaker as more data is collected. 
A separate study does not find any correlation [2] . In ad- 
dition, a Greizen-Zatsepin-Kuzmin mechanism (GZK), 
described as a photo disintegration of the nucleus on cos- 
mic microwave background, see [3], limits the UHECRs 
travel distance, effectively excluding the distant sources. 
A flux suppression that could be interpreted as the GZK 
cutoff was independently seen by two large cosmic ray 
observatories [4, 5]. Also, current theoretical models, 
with the exception of the most exotic "top down" mod- 
els, fail to predict mechanisms capable to accelerate the 
cosmic rays above 10 20 eV, yet the cosmic rays above 
that energy were observed [6]. The controversy in the 
chemical composition of the UHECRs has persisted for a 
long time [7, 8]. Despite so many unanswered questions, 
the UHECRs remain to be the only source of elemen- 
tary particles with energies above the reach of modern 
accelerators and allow us to study particle interactions 
at these energies [9]. In addition, the UHECR research 
can be considered as "non-electromagnetic" astronomy, 



complementing research in different electromagnetic fre- 
quency bands. 

The cosmic ray spectrum falls rapidly with energy, and 
the very low cosmic ray flux makes direct observations 
at the energies above 10 14 eV impractical. In order to 
increase the detector volume, the earth's atmosphere can 
be used as a giant calorimeter. Cascades of secondary 
particles in the atmosphere (extensive air showers) can 
be detected by air fluorescence or Cherenkov radiation in 
UV and by particle counters on the ground. The modern 
hybrid cosmic ray telescopes utilize two or all three 
of the above mentioned techniques simultaneously [10, 
11]. The air fluorescence technique, which is capable of 
observing a large portion of an extensive air shower in the 
atmosphere, provides very good primary particle energy 
resolution combined with a precise pointing to the arrival 
direction. Detailed information about the shower profile, 
the number of charged particles as a function of the slant 
depth in the atmosphere, and the depth of the shower 
maximum (X max ) in particular, can be used for chemical 
composition and particle cross-section study [7, 8, 9]. 
The UV observations, however, can only be conducted 
during moonless nights, which limits the air fluorescence 
and the air Cherenkov detectors duty cycle to not more 
than 10%. The ground particle counters can be operated 
continuously, increasing the data statistics, but the energy 
and the shower profile reconstruction for this method 
relies on extensive Monte Carlo (MC) simulations, and 
is not yet as accurate as the air fluorescence technique. 



Detection of the extensive air showers caused by 
UHECRs by radio emission was first used in the 1960s 
[12]. The success however was somewhat limited due 
to a lack of understanding of the radio emission mecha- 
nism. The later efforts were concentrated more on the air 
fluorescence, air Cherenkov and the ground counter ar- 
ray methods, but the radio observation technique is very 
appealing with such an advantage as 100% duty cycle 
combined with lower deployment and operational costs. 
The advancements in modern computing for the exten- 
sive MC simulations led to a better understanding of 
the radio emission mechanism and the new ground and 
balloon-borne experiments are in operation or planned 
[13, 14, 15, 16, 17]. 

Similar to a ground counter array, the energy and the 
X max reconstruction from radio data relies on the MC 
simulations, which, in turn, require a complete under- 
standing of the emission from a moving charge in a dense 
medium. A "full" MC approach tracks radio emission 
from each individual particle and takes into account co- 
herence effects of the air shower for a ground observer 
[18, 19]. The results of the MC simulations derived us- 
ing different approaches start to agree with each other, 
but only in the far zone and away from a low frequency 
limit [20]. In this work, we compare the different theo- 
retical approaches to the problem of the emission from a 
charged particle in a dense medium to study their validity 
and limitations. 



taken into account, assuming that the particle does 
not emit while moving uniformly in vacuum [26]. 
The medium effects are taken into account at the 
end by adding the index of refraction. 

To derive both approaches, let's consider a charged par- 
ticle moving along the track %(x' s ), see the diagram on 
Figure 1, where apostrophe indicates a delayed time or 
coordinate. 
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FIGURE 1. Radio emission from a source in a dense 
medium. 



2.1. Summing over the straight segments of 
the track. 



2. EMISSION FROM A MOVING 
CHARGE IN A DENSE MEDIUM 

Despite appearing to be a "classical" problem, the elec- 
tromagnetic emission from a moving charge in a dense 
medium does not have a complete analytical solution 
yet. The problem was first formulated by Frank and 
Tamm in 1937 [21], but was only solved for some special 
cases such as the far zone approximation, infinite particle 
tracks, slowly changing properties of the medium, a par- 
ticle passing through a thin plate and a few others [22]. 
For an infinite track in vacuum the classical solution can 
be derived from Maxwell's equations with no difficulty 
[23, 24]. In case of a finite track in a dense medium with 
a constant index of refraction n > 1 the track can be split 
into small straight segments and then two different ap- 
proaches can be used: 

• only the emission from the straight part of a seg- 
ment, where the particle moves uniformly in a dense 
medium, is taken into account. The emission from 
the ends of the track is ignored assuming that the 
emission time is very short [25]. 

• only the emission from the ends of a segment, 
where the particle accelerates or decelerates, is 



Let's first express electric and magnetic fields using 
electric and magnetic potentials: 

B = VxA (2) 
where the delayed potentials at time t = t— r ^ Xs ^ s > are: 

^ f ) = i^/-r^4 dV ' (3) 

4n£o J r(x s ,x s ) 

A(x s j) = -^jj^dV (4) 
47T£oc 2 J r(x s ,x s ) 

and p(x s ,t) and j(xj ,t) are the charge and current densi- 
ties respectively, and x s and x s are the particle coordinate 
in the observer and the source reference frames. With the 
Fourier components of the fields defined as: 

E(x s ,t)=Re J E (0 (x s )e~ im dco (5) 
B(x s ,t) = Re [ B (D (x s )e- im dco (6) 



the electric and magnetic field spectral components can 
be written as: 
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where k = ®f, r is a unit vector along vector r from 
the observer to the particle coordinate, r = \f\ and the 
spectral components of the charge and current densities 
are: 
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The common approach is to drop the static 1/r 2 term 
when differentiating eq. 7, but we will have to impose 
an even stricter condition later. After differentiating eq. 7 
we obtain: 
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In the far zone, where kr » 1, E _L B, the spectral 
component of the electric field can be written as: 
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or from eq. 11: 
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In the Fraunhofer zone R » \x(t')\ during the whole 
observation time, where R = |/?| and R is a vector from 
the observer to the source coordinate origin. Also r{t') ~ 
R — R ■ %{t'), where R is a unit vector along R, thus, we 
can replace r by R and r by R in eq. 12. For a charge q 
moving along the trajectory 2(0' me current density can 
be expressed through the Dirac delta function: 
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Substituting eq. 9 and eq. 13 into eq. 12 and taking into 
account that j8 = and k = n®r one can obtain 

the expression for the spectral component of the electric 
field: 
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where /3± = fx/5 x r is the perpendicular component 

of the charge velocity directed along R. 

For the particle acceleration up to time to, see the dia- 
gram on Figure 2, the spectral component of the electric 
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FIGURE 2. Particle velocity while moving along a track 
segment. 



field is: 
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The integral is eq. 15 can be approximated by f3sinOto, 
where 6 is the angle between j3 and R. Since j3 and sinO 
cannot exceed 1, the whole integral can not exceed £r> 
If << A^, we can safely neglect the emission from 
that end of the track segment. Similar conclusions can 
be made for the deceleration end of the track segment. 
Integrating the remaining part of the expression 14 and 
summing over all track segments we obtain: 
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where m is the segment index. 



(16) 



2.2. Summing over the emission from the 
ends of the track. 

In this approach we account for radiation from the 
ends of the track %(x' s ), see the diagram on Figure 1, 
where the charged particle accelerates and decelerates, 
and ignore the radiation from a straight portion of the 
track where particle moves without an acceleration, as- 
suming the particle moves in vacuum first and adding 
the index of refraction later. We start with the equation 



for the electric field [23]: 
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where all values are taken at a delayed time t' . In the 
Fraunhofer zone R « r, see Figure 1 , and magnetic and 
electric fields are perpendicular: B(t,x s ) = ^E(t,x s ) xr. 

Also note that = 1 — rft For the spectral component 
of the electric field we obtain: 
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we finally obtain for the electric field: 
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Solving eq. 20 for particle birth or acceleration end of the 
track segment, where t < to, yields: 
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and for the deceleration end, where t > to + At, yields: 
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We can neglect how the particle velocity change during 
acceleration or deceleration if coSt « I. Finally, adding 
the emission from both ends of the track segment, ac- 
counting for the index of refraction and summing over all 
segments of the track we arrive to eq. 16 again. Note that 
the index of refraction was added to the solution for the 
spectral component of the electric field at the very end. 
Strictly speaking this is not correct because we started 
with the full spectrum equation for the electric field 17. 
In a narrow frequency band however, where the medium 
dispersion can be neglected, such an approximation is 
valid. 



3. COHERENCE WAVELENGTH AND 
COHERENCE TIME. 

We split the particle trajectory x( x s) into Fresnel zones 
so that the phase difference between nearby zones 
is equal to 71, see the diagram in Figure 3, so that 




FIGURE 3. Particle trajectory split into Fresnel zones. 

ja*{l-k 2 &)+in = e im(l-RSi)^ wne re R x and R 2 is the 
distance from the observer to the corresponding Fres- 
nel zone, and ft and ft are corresponding velocities. 
In the Fraunthofer limit R\ w R 2 « R, ft « ft « j8 and 
At = t 2 - h • Then At(o(\ Rfi) = n or 
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note that 
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and Rf3 = ft cos 6, where 6 is the 



observation angle. The left hand side in the eq. 23 has 
the dimension of length, which is called the coherence 
length: 1 
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The integration along the whole track can be replaced by 
the integration over one Fresnel zone. For example, using 
the Fresnel zone closest to the observer for the emission 
by the particle coming from infinity and stopping at time 
to'. 
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The choice of the Fresnel zone is arbitrary. Although 
the particle is emitting from the whole track, it appears 
for the observer as if all the emission comes from the 
track's end. The emission from the nearby Fresnel zones 
is canceled out due to the phase being different by 7t, 



Since eq. 24 is not valid at the Cherenkov angle, 6 C , another expres- 



sion should be used: L co h - 
by I.Tamm in 1939 [27]. 



, Related definitions were formulated 



with exception for the closest, or any other arbitrary 
chosen Fresnel zone, where an asymmetry is imposed. 
One practical application of this is the reduction of the 
computing time for the MC simulations of the radio 
emission from air showers. 



4. CONCLUSIONS. 

Both approaches presented above are mathematically 
equivalent if both derived under kf » I condition. In a 
low frequency limit as well as close to the Cherenkov an- 
gle, the second approach can go to infinity, while the first 
give a finite answer, but both approaches, strictly speak- 
ing, are not valid. One also should be careful about sub- 
tracting very close numbers doing MC simulation on a 
computer when using the second approach. As noted be- 
fore, adding the effects of the medium into the second 
approach is not straightforward, but a simple addition of 
the index of refraction is valid as an approximation for 
some special cases. Another interpretation of the prob- 
lem is proposed in [28]. The MC simulation codes de- 
scribed in [19] and [18] are based on the first and second 
approaches respectively. The results of these MC simu- 
lations have never been validated by an accelerator ex- 
periment. The full solution without the far zone limita- 
tion can be obtained, but the resulting integral contains 
a modified Bessel function of the second kind and can 
not be solved analytically. Finally, it should be noted 
that we should interpret the emission in the presence of 
medium as emission from all moving particles, includ- 
ing the charges in the medium. We did not account for 
the particle interaction with the charges in the medium 
in both approaches. This is not a problem for the MC 
simulations on the computer. 
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